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Abstract 

(N ^ _ . . + 

^ I The process of e^e annihilation into a tt^vt pair with radiation of a photon is con- 

' sidered. The amphtude of the reaction e^e" — > 7r^7r~7 consists of the model independent 

initial-state radiation (ISR) and model dependent final-state radiation (FSR) . The general 
structure of the FSR tensor is constructed from Lorentz covariance, gauge invariance and 
^ discrete symmetries in terms of the three invariant functions. To calculate these func- 

. tions we apply Chiral Perturbation Theory (ChPT) with vector and axial-vector mesons. 

\ The contribution of e'^e~ — > 7r"*"7r~7 process to the muon anomalous magnetic moment is 

Q^' evaluated, and results are compared with the dominant contribution in the framework of 



Oh! 



a hybrid model, consisting of VMD and point-like scalar electrodynamics. The developed 
(D ', approach allows us also to calculate the vr+vr" charge asymmetry. 



JXl ■ PACS: 12.20.-m; 12.39.Fe; 13.40.-f; 13.66.Bc 

1 Introduction 

The cross section of electron-positron annihilation into hadrons, e^e~ — > hadrons, is crucial 
for evaluation of the hadronic contribution to anomalous magnetic moment (AMM) of the 
muon a|j"'^, and is at present one of the main sources of theoretical uncertainty in the value of 
AMM p. In order to resolve the existing deviation of the experimental and Standard Model 
prediction values of AMM, the corresponding hadronic contribution is needed with very high 
accuracy, better than 1%. This is especially important in view of a new E969 experiment, which 
is expected to measure AMM about three times more accurately the updated £'821 experiment, 
in which the data are expected to have AMM about three times more accurate than now 
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The hadronic contribution to AMM cannot be reliably calculated in the framework of per- 
turbative QCD, because low-energy region dominates, and one usually resorts to dispersion 
relations, where the experimental total cross section is the input. Experimentally, the energy 
region from threshold to the collider beam energy is explored at the ^-factory DA$NE (Fras- 
cati) j3] and S-factories BaBar (SLAG) and Belle (KEK) ^ |H] using the method of radiative 
return [HI [7| |H]. In spite of the loss in the luminosity, this method potentially may have ad- 
vantages in systematics over the more traditional direct scanning measurements performed at 
different CM energies, such as experiments on VEPP-2M (Novosibirsk) 9] and BES (Bejing) 

cni. 

The radiative return method relies on a factorization of the e'^e hadrons + 7 cross 

section in the product of the hadronic cross section a{e^e~ hadrons) taken at a reduced 
CM energy and a model-independent radiation function known from Quantum Electrodynamics 
(QED) [HI El El- This factorization is valid only for photon radiation from the initial leptons 
(initial-state radiation (ISR)). The additional contribution from photon radiation off the final 
hadrons (final-state radiation (FSR)) is model dependent and becomes a background in the 
radiative return scanning measurements. That is why the problem of the separation of ISR and 
FSR has become quite important. 

Different methods have been suggested to separate ISR and FSR contributions for the 
dominant hadronic channel at low energies, mainly the pion pair production process 

e"(pi) + e+(p2) ^ vr+(p+) + 7r-(p_) + 7(fc). (1) 

One of them is to choose a kinematic set up, where the photon is radiated along the momenta 
of the leptons (DA$NE setup, [3 E] and references therein). In these conditions the FSR 
contribution is suppressed. If the FSR background can be rehably calculated in some theo- 
retical model, then it can be subtracted from the experimental cross section of process (Q) or 
incorporated in the Monte Carlo event generator used in the analysis. Finally, the theoretical 
predictions for FSR can be tested by studying the C-odd interference of ISR and FSR jSl E] • 

Another, and even more important reason why one should know the FSR cross section is the 
fact that the next-to-leading order hadronic contribution a^^'^''* to AMM, where an additional 
photon is attached to hadrons, is of the order of the present experimental accuracy. 

The FSR cross section in process (P) has been calculated El in the framework of 
scalar QED (sQED), in which the pions are treated as point-like particles, and the resulting 
amplitude is multiplied by the pion electromagnetic form factor F^(s) evaluated in the Vector 
Meson Dominance (VMD) model (s is the invariant e+e~ energy squared) to account for the 
pion structure. In this model the contribution of the channel 7r"'"7r~7 to AMM is estimated as 
^hada ^ 4 3 ^ ^Q-10 Although sQED in some cases works well [3 El, it is clear that sQED 
is a simplified model of a complicated process, which may include excitation of resonances, loop 
contributions, etc. In view of the above mentioned requirements for the accuracy of theoretical 
predictions, further studies of the FSR contribution are necessary. 

In this paper we consider the e"'"e^ — 7r"'"7r~7 reaction in detail, focusing on FSR. Firstly, we 
specify the model-independent structure of the FSR amplitude, based on Lorentz covariance, 
gauge invariance and discrete symmetries. Taking this structure into account we rewrite the 
FSR contribution, as well as the interference of ISR and FSR, in terms of the three scalar func- 
tions fi depending on three kinematical invariants. Secondly, the model-dependent functions fi 
are obtained in the framework of Chiral Perturbation Theory (ChPT) with vector p(770) and 
axial- vector ai(1260) mesons [14j. In this way fi are expressed through the several constants 
entering the ChPT Lagrangian. 
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Figure 1: Diagrams describing e+e tt+tt 7 process 



For experimental conditions, in which the cross section integrated over the full phase space 
of the two pions is required, this integration is carried out analytically. We obtain a general 
result for the cross section da/dq'^dcos9 in terms of the two scalar functions hi^2{,(f, s) (g^ is 
the invariant mass of the 7r+7r~ pair, 9 is the angle between photon and electron momenta). 

We further study the interference of ISR and FSR by calculating the 7r"'"7r~ charge asymme- 
try. Finally, the contribution for '7r"'"7r~7 channel to a^^'^'"' is evaluated and results are compared 
with sQED predictions. 

The paper is organized as follows. In Sect. |21 the general form of the amplitude is considered. 
The squared and averaged amplitudes for ISR and FSR, as well as the interference part are 
analytically calculated and the structure of the cross section dcr/dg^dcos6' is studied. The 
invariant functions fi in framework of ChPT are derived in Sect. El Results of calculations and 
discussion are presented in Sect. El In Sect. Owe draw conclusions. In Appendix A we discuss 
symmetries of the FSR tensor and its model-independent structure. Appendix B contains the 
ChPT Lagrangian and explicit expression for the FSR tensor. In Appendix C the Feynman 
rules needed for evaluation of the 7* — * 77?+ vr" amplitude are specified. Appendix D collects 
expressions for scalar functions hi 2{q'^,s) in ChPT. Some subtle aspects of the kinematics at 
low values of are presented in Appendix E. In Appendix F the contribution to the FSR 
tensor and charge asymmetry from the "anomalous" p ^ 717 process is calculated. 



2 Formalism for e+e tt+tt 7 reaction 

Reaction (1) is described by diagrams depicted in Fig. 1. To analyze it we introduce 4- 
momenta Q = Pi + P2, Q = P+ + P- and I = p+ — p-. The amplitude of the reaction depends 
on five independent Lorentz scalars, which can be chosen as follows: 

s = = 2pi- p2, h = {pi - ky -ml = -2pi ■ k, t2 = {p2 - kf - ml = -2p2 ■ k, 

Ml = l-pi, U2 = 1-P2, (2) 

where we neglected the electron mass (mg) in the expression for s. For further reference note 
that other invariants are related to those in Eqs.(j2I), f*^^ instance, = s + ti + t2, = 
Ami - s - ti - t2, ? ■ ^ = 0, k ■ Q = k ■ q = -^{ti + t2), k ■ I = Q ■ I = ui + U2. 

The amplitude of process (|H) is the sum M = Mjsr + MpsR, where Mjsr (Mpsr) describes 
ISR (FSR) process 

Mjsr = ^L^'ell^F^iq^), MpsR = -J.M^'el, (3) 
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where the lepton currents are given by 



J., 



e'us^{-p2)[Y— ^ -Y + 1 -i]usApi 



(4) 
(5) 



FT^{q^) is the pion electromagnetic (EM) form factor, e* is the photon polarization vector and 
the tensor Mp^ describes the photon radiation from the final state. This tensor is considered 
in detail in Sect. El and Appendices A and B. In Eqs. (jl)) and © Us^{pi) and 'Us2(— P2) are the 
electron and positron spinors with normalization: Us'{p)Us{p) = ~Us'{—p)Us{—p) = 2me5ss'- 

The invariant amplitude squared, averaged over initial lepton polarizations and summed 
over the photon polarizations^ can be written as 



IMI 



\Misr\'' + WpsRi" + 2Re(M,5i?M 



FSRJ- 



(6) 



The expressions for IMisr]"^, IMfsrI"^ and the interference part Re{MjsRMpgj^) are given in 

Sects, lino and E31 

The differential cross section for process is written in the following form 



2E, 2E_ 2uj 



IMP 



(7) 



where p+ = (£■+, 



{E_, p_) and k 



|k|,k). 



2.1 Initial- state radiation 

Let us consider first the ISR contribution shown in Fig. la. The amplitude squared can be 
written as 



\MjsR\' = -^jF^iq')\'L^;jFr, 



(8) 



+ 



tit2 



2' 



''1 



''2 



6m 



''1 



\P2^,P2u 



tlt2 



6m; 



where pi^ = pi^ - q^{pi ■ q)/q^ and similarly for pa^, and g^^ = g^,^ - q^,qu/q^. 

If one integrates over the full, unrestricted phase space of the final pions, the hadronic tensor 
can be integrated in the invariant form (see [TT] ) 



j 5^{pi + p2 - p+ - p- - k) l^l^ 



d^p+d^j9_ q 



which leads to the following cross section ( "F" stands for "Full" ) 



2 i(^{<f)9t,u, 



(9) 



da 



(F) 
ISR 



dcudl] 27r2s 



a{q^ 



tit 



1^2 



ti 



(10) 



iWe use Epo/ar. = -9pa 
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where criq^) = — is the total cross section for e~^e^ tt+tt^ (), C = (1 — 4m?./g^)^/^ 

and a = e^/Air ^ 1/137 is the fine-structure constant, dQ = dcos^^d0 and 6 (0) is the polar 
(azimuthal) angle for the emitted photon. Note that the cross section does not depend on the 
azimuthal angle (p. 

Integration over the photon angles leads to the result 

!!^ = .fe^)i(^±4(L-l,, L = ta4. (11) 

dq^ ns^{s — q'') mj 

In the last equation we used to = (s — g^)/2ys and dg^ = —2^/sduJ. Note that Eq. (fTTj) holds 
for the full angular phase space of photon; another case of the restricted angular phase space 
of photon will be studied elsewhere. 

In some cases, due to experimental conditions, the entire phase space of the pion is not 
available. Then it is not possible to use Eqs. ©-(HH)- In this situation one has to contract first 
the hadron and lepton tensors and then carry out phase-space integration. From Eq. (jH)) we 
obtain the ISR contribution 



\Misr\' = ^\F^{q')\''R, (12) 



m. 
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R = -^F + 



Xl + xl~ Xiiq'^ - h) - X2{q^ - ti) 



'^ml (xi _ \ _ 2m| / X2 _ ^ 



tl \q^ ) tl \q^ 

where Xi,2 = 2pi,2 ■ V- = |^i,2 ~ ^-^ ~ ^i,2- Then the ISR cross section takes the form jT2I ("R" 
stands for "Restricted") 

da\% = ]M^afkd^\p^]dE^ - ^ - - E_), (13) 

where we introduced the electron (positron) energy E = y/s/2 in the CM frame. Using the 
relation 

np^ldE^dc^ r IP+pdc^ 

J E_ ^ y {2E -Lj)\p+\+uE+cose^+ ^ ' 

we obtain the corresponding ISR cross section 

_ Mq') « ^ }^^^^^ 7d^ T R\p^\^dc^ ^^^^ 



dg2 4E^( 27: E J J 271 



If < 2E{E — mT,)/{2E — m^) (which, for example, corresponds to values q^ > 0.16 GeV^ at 
^/s = 1 GeV) then we obtain 

_2E{E- uj) {2E -uj)-uj cos 0^+ ^AE^{E - - ml [{2E - ujf - uj'^ cos^ 6^+] 

^+ ~ (2E-CU)2-CU2C0S2 9^+ • ^^^^ 

Here c+ = cos 6'+, 9+ (v?+) is the polar (azimuthal) angle of the positively charged pion (we 
take OZ axis along the vector pi), and cos^-y+ = kp+/|k||p+|. In this energy region the angle 
6'^+ can take arbitrary values (for details see Appendix E). Other notation can be found in 

Ref. ng. 
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2.2 Final-state radiation 



The process of photon radiation from the final state is shown in Fig. lb, where the dark 
rhomb denotes the set of the contributing diagrams. The covariant decomposition of the FSR 
tensor M'^^ can be obtained from Lorentz and discrete symmetries (Appendix A). This tensor 
involves three gauge- invariant structures and invariant functions /, {i = 1, 2, 3). The explicit 
form of the functions in the framework of ChPT is discussed in Sect. El In terms of fi we 
obtain 



FSR\ 



%\FJs)\' 



an|/ip + 2ai2Re(/i/*) + 0221/2! 



+ 2a23i?e(/2/3*) + asal/sl' + 2arsRe{hf; 



(17) 



an 

^22 



1 

-S{tl + tl), 033 = - — [tit2p + 2(mi + U2){u2ti + ^1^2)], 



sl\h + t2f + Al' 



ui\s^ + s{ti + t2) + ta) + U2 + s{ti + ts) + ti) 



+2MiM2[s^ + s(ti +^2) -M2] 



+ s{u\ + u\){ui + U2) 



ai2 



ai3 

^23 



sl\ti + ^2) + 4mi"(s' + st2 + ta) + 4m2 (s' + sti + t{ 



+AuiU2[2s^ + s(ti + 12) - 2tit2 



-[{Ui + ■U2)(sti + St2 + ^1^2) - Ultl - U2tl] 



-l'^{uitl + U2t\ - {ui + ^2)^1^2) - 2s{ui + ^2)^ 

+2{ui + 'U2)[m1^^2 + M2^il - UiU2{tl + ^2)] 



In order to obtain the cross section dapsR we have to substitute Eq. (fT7|) in Eq. (|7|) and 
integrate over the phase space of the final particles. 

In the case of the full phase space of pions, the integration can be simplified using the 
method suggested in Ref. [TT]. In this case the squared matrix element \Mfsr\'^ in Eq- © can 
be integrated in the invariant form 



\Mfsr? 



E^E+ 

Taking into account the conditions W^^Q^ = W^^Q^ = 0, one can write W^^ as 



W^"" = hig 



where hi^2 are functions of and . 



^(/ii + g2/,2)A;^r + h2Q^Q' - ^i + ^'^2 (fcMg^ + k'^Q^) 



k-Q 



(19) 



6 



In the framework of sQED /ii 2 were found in Ref. ^T] (see also Appendix D). Using Eq. (jl9j) 
we obtain the following equations for /ii 2 in terms of the functions fi for any model of FSR : 

h2{k-qf = j^^^S\Q-p^-P--mk.lf\h\^ + {k.Qf\h\' 

-2{k ■ l){k ■ Q)Re{hfm{k ■ Qfl^ + {k ■ lf{Q^ - 2k ■ Q)], (20) 

2h-h2Q' = j'^^^5\Q-p^-p^-k){2{k.Qf\h\' 

+ [l\k ■ + I'ik ■ lf{Q^ -2k-Q)+ 2{k ■ OII/2I' 
+2[{k ■ ifQ^ + l\k ■ Qf]Re{hf*) - AQ\k ■ Q){k ■ l)Re{hf^) 
-4Q\k ■ ifRe^hfl) + [(^ ■ O'Q' + 2(fc ■ lfQ\k ■ Q) 
-fQ^k.Qfm'} (21) 

with Q"^ = s. In Appendix D the explicit form of hi^2 is presented in the framework of ChPT. 
Then the cross section dapg^/dujdcosO takes the form 

da^fsR _ '^^^ r, tit2 



Integrating this equation over the polar angle of the emitted photon we find 

If we deal with the restricted phase space we can use the same arguments which led to 
Eq. ()15|1 in Sect. 12.11 Then the cross section can be written as 

dg2 32s(27r)3v^(2E-cu)|p+|+cuE+c^+ ac+acosa, 



where \Mfsr\'^ is determined in Eq. (fT7j) . 



2.3 Interference 



The interference part of the squared invariant amplitude |M|2 is written in terms of the 
invariant functions fi and the pion form factor as follows 



Re{MisRM*psj^) 



where the coefficients Aj are 



AiRe(F.(g2)/*) + A2Re(F.(g2)/*) + A.ReiF^iq') f; 



(25) 



A O + ih-t2)s + ht2 + + {t2-h)s + tit2 , + 

Ai = -2Ui{ + — )+2M2( + — ), 

ti 12 12 ti 



A2 = -4!^i!±^ 

tit2 



{s{tl - t2) - t2{h + t2))ui - {S{t2 - tl) - tl(tl + t2))u2 



tlt2 



Ui[2s'{ti - t2) + Siti + t2){ti - 3^2) - t2{ti - t2y 



- U2[2s\t2 - h) + S(tl + t2){t2 - 3tl) - h{t2 - t 



I) J h 



A, = -2.(ti-t2)/^-4.( V(2 + ^)-V(2 + ^) + %-^nin2). (26) 

ti 12 tlt2 
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We would like to mention that the cross sections dajsR and dapsR are symmetric under the 
interchange of tt"*" and vr" momenta, while the interference term dajpR is antisymmetric: 

daisR ip+ , p- ) = daisR {P-,P+), dapsR ip+ ,P-) = dapsR {P- , P+ ) 
daiFR{p+,P-) = -daiFRiP-,p+). (27) 

Therefore dajpR integrated over the symmetric phase space of the pions (for example, for the 
full unrestricted phase space) is equal to zero. Other implications of Eqs. (|77j) are considered 
in the next sections. 

For the restricted phase space of pions we have a result analogous to 

d^ = MMjsrM^sr) - JP±l ^dcdcos^. (28) 

dg2 16s(27r)3 ^ {2E - uj)\p+\ + u;E+c^+ 27i + ^ ^ 

3 Final-state radiation in the framework of ChPT 

Based on results of Appendix A we can write the FSR tensor Mj"^ in the form 

Mr = /irr + /2rr + fsrr- (29) 

In framework of ChPT with vector and axial-vector mesons ^1] (see details in Appendix B) the 
process 7* 7r+7r^7 is described at tree level by the diagrams shown in Fig. 2. Using results 
from Appendix B we find the invariant functions fi = fi{Q'^, k ■ Q,k ■ I) 

U = ff^"" + A/„ (30) 

where f^^^^ correspond to sQED 

fsQED _ 2k-QF^{Q^) 
Jl 



{k-QY-{k-l) 



2 ■ 



fsQED _ 2F^(Q^) fsQED _ f) /oiN 

J'2 ~ 7Z 7yv2 7Z a2' ~ ^' ) 



{k-QY-{k-lY 
and A/j are additional contributions 

- 2FvGv f I 1 



A/i 



F 



/2 \ml ml - Q2 

2 + (^-0' , {Q^ + k-Q)[Aml-{Q^ + l' + 2k-Q)] 
D{l)D{-l) 8D{l)D{-l) 



(32) 



_ _ n ^ml-iQ' + l' + 2k-Q) 

- nml2D{l)D{-iy ^^^^ 

Here D{1) = ml - {Q^ + + 2k ■ Q + Ak ■ I) /4. The functions Afi and f-^^^ obey the same 
symmetry relations as given by Eq. ()A.5|) for functions fi. The EM form factor in Eqs. (jHT|) for 
the on-shell pion follows from Eq. ()B.4|) : 
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Figure 2: Diagrams for FSR in the framework of ChPT. Dashed hnes depict pion, wavy hnes 
- photon, double-dashed hnes - p"— meson, and dotted hnes - Oi— meson. The hatched circles 
denote the irreducible 77777 vertex 
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meson 


m (GeV) 


Gv (GeV) 


Fv (GeV) 


Fa (GeV) 


P 


0.775 


0.066 


0.14 




ai 


1.23 






0.122 



Table 1: Masses and coupling parameters of vector and axial- vector mesons 



To account for the finite width of the vector meson one can substitute in Eqs. ()32j) and ()35p 
ml-Q'^ml-Q^- ^m,T,{Q% T ,{Q^) = ^^^^(1 - ^-^?''0{Q' ' ^O. (36) 

where rp(Q^) is the energy-dependent width for the p ^ -n-n decay. In a similar way one can 
include in Eqs. (jn^ -()34j ) the width of the decay 37r. The analytical form of Y can 
be taken from, e.g.^ Ref. 

Using the form of /i.2,3 we can find the functions h\^2 from Sect. 12.21 appearing in the FSR 
cross section. The expressions are rather lengthy and are listed in Appendix D. 

We would like to mention that the Compton 777 77? scattering amp litude in the frame- 
work of ChPT 14j was calculated in Ref. [TE]. Having compared Eqs. (37)-(40) of Ref. [12] 
with Eqs. ((221)- (EH) of this paper, we have concluded that the p-meson contributions are equal 
whereas the ai-meson contributions are different. 

In addition to the even-intrinsic-parity contributions considered above there is an odd- 
intrinsic-parity part. The corresponding Lagrangian f7[ QHl describes processes which do not 
conserve intrinsic parity, such as p — 7r7. The contribution of the two-step mechanism 7* — >■ 
p^Ti^ — > 7r+7r~7 to the FSR tensor is evaluated in Appendix F. 



4 Results of calculation 

Table 1 lists the parameters of the model. The couplings Fy^F^ are determined from the 
experimental decay widths [19 : r(p'^ e^e~) = 6.85 ± 0.11 keV and r(ai 717) = 640 ± 246 
keV, while Gy is fixed from the width r(p —>■ tttt) = 150.7 ± 2.9 MeV (we neglect the chiral 
corrections here). The pion weak-decay constant is f-j^ = 92.4 MeV. 



4.1 Charge asymmetry 

We will illustrate the results obtained in the previous sections by considering the charge 
asymmetry jTHI 

^- Nie^) + Nie^) ^^^^ 

for "collinear" kinematics in which the hard photon is radiated inside a narrow cone with the 
opening angle 26*0, 9o -C 1, along the direction of initial electron. We choose 9o = 7°. In these 
conditions the asymmetry takes the form 



A 



da 



{R} 

IFR 



dg^dc-i 



da 



ISR 



dg^dc_| 



-1 



(38) 



where we neglected the FSR contribution compared to the ISR one in the denominator. 
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Figure 3: Charge asymmetry as a function of pion polar angle at fixed for s = 1 GeV^. Here 
the solid line corresponds to sQED, the dashed line - the full result in ChPT. 
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Figure 4: Charge asymmetry as a function of at fixed pion polar angle for s = 1 GeV^ 
Notation for the curves is the same as in Fig. 3. 



The ISR cross section for collinear kinematics was obtained in Ref. ^21- We use Eqs. (26)- 
(30) of Ref. (t2J, which were derived in the quasi-real-electron approximation. It is convenient 
to rewrite these results as follows 

d4|L^ .a^^|F,tf)P .^ . (39) 

P{z, Lo) = — ^Lo ^, Lo = In ■ 



s{s — q^) s — q^' Ami 



e 

2 ^2 „2 



Vl z\{\^z)K-{\-z)c^'U rrXi_xi_< ^ q_ 
i K[{l + zY-{l-zYclf' s s2 g2 ' ^ 5' 



In order to obtain dcrj^^^/dg^dc+ we should integrate the right-hand side of Eq. (PHj) over 
<y9+ and 9. Since the right-hand side has no singularity at the point 6* = the integration over 
(p+ and 9 can be easily done numerically. 

In Figs. 3 and 4 we show the asymmetry dependence on pion polar angle (at fixed q^) and 
on q^ (at fixed pion polar angle). It follows from the calculations that the asymmetry changes 
sign at about q^ = 0.5 GeV^ (see Fig. 4). At all pion angles the difference between sQED and 
ChPT shows up only at small values of q^, i.e at the high photon energies. 

The sQED description is adequate for soft photon emission. It follows from Eqs. (jHTjl - ljH^ 
that at small photon energies, f^'^^^ ~ ff^^^ ~ whereas A/j which are responsible 
for the deviation from sQED behave rather as constants. Only at large photon energies the 
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Figure 5: Differential contribution to a/j " ''^ from e+e tt+tt 7, where 7 is the hard photon 
with energy uj > E^ut (left panel). Integrated contribution to function of Ecut for 

Smax = 1-5 GeV^ (right panel). The solid (dashed) line corresponds to the sQED (full) result. 



contribution from the intermediate Oi-meson (the last two diagrams in Fig. 2) can be sizable, 
because the denominators D{l)D{—l) in Eqs. approach the ai-meson pole with the 

photon energy increasing. 

For high value of the difference between predictions of sQED and full calculation in ChPT 
is small. For example, at = 0.8 GeV^ and = 0.6 GeV^ it is less than 1% (the dashed and 
solid lines in Figs. 3 and 4 almost coincide). Taking into account that the asymmetry itself is 
less than 10~^, the experimental observation of such deviations in the energy region > 0.6 
GeV^ is problematic. 

Additional contribution coming from the process 7* — >■ p^ir^ —>■ 7r+7r^7 turns out very small 
(see Appendix F) and does not change the above conclusion. 

In order to test the calculation we can check that the asymmetry, integrated over the 
symmetric phase space of the pions, vanishes. Since no restriction has been imposed on the vr" 
polar angle we impose no restriction on the vr"*" polar angle, i.e. choose 0° < 6^,6^ < 180°. 
Indeed, the integrated asymmetry is equal to zero. 



4.2 Contribution from e+e tt+tt 7 to AMM of the muon 

Now we apply the previous results to the calculation of a^'^'^''^ which arises from the 7r+7r^7 
channel. We should mention that only the radiation of the hard photon with the energy 
^ > Ecut is taken into account. 

According to [201 a^^'^'^ can be written in terms of the dispersion integral 



K{s) = C —^iL^L-da;, 
Jq + [I - x)s/mf^ 



where is the muon mass. 



To obtain the cross section a'^ '"^(s) we have to integrate Eq. (j23|) over q^ from 4m^ 
to q^ax- The value of q^^x determined from the equahty q^^x = ^ — 1^E(,ut- From the 
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0,2 0,3 
E„.(GeV) 



0,4 



Figure 6: Integrated contribution to a^'"^'^'^'' as a function of Ecut for Smax 
for the curves is the same as in Fig. 3. 



1.5 GeV^. Notation 



condition g^^^ > 4m^, the lower hmit of the integration in Eq. pUj) is found to be Smin = 
{Ecut + ^/4:m'^ + Ecut^- The upper hmit in pUj) is replaced by a finite Smax- The value of Smax 
is chosen s^ax = 1-5 GeV^, which is about of 0(m^J, the upper limit of the applicability of 
ChPT with p and Oi mesons. The dependence of a'^f''^''^ on energy E^ut is shown in Fig. 5. 

As follows from our calculations the additional contributions to a^"''^''^ stemming from ChPT 
are very small compared to the sQED result. This is in line with the conclusion of Ref. |21j . 
Even for relatively large cut-off energy Ecut = 200 MeV the full result in ChPT differs from the 
sQED result only by 3.5%. For that reason the solid and dashed lines in Fig. 5 almost coincide. 
At the same time with increasing photon energy sQED begins to loose its predictive power. 
This is demonstrated in Fig. 6. In this region of energies the contribution from the ai-meson 
is considerable and has to be taken into account. For example, at a photon energy about 500 
MeV the deviation from sQED reaches 60%. However such a deviation (which is of the order 
of 10~^^) is beyond the accuracy of the present measurements of the muon AMM. 



5 Conclusions 

In this article the FSR of a hard photon in the e'^e" — > 7r+7r^7 reaction has been considered 
in framework of ChPT with vector p and axial-vector ai mesons. The respective Lagrangian 
generates effective chiral terms and, as substantiated in jHj, is adequate for description 

of processes at energies up to about 1 GeV. 

Our consideration of FSR is motivated by the necessity to study the model dependence 
of the hadronic contribution a^"-'^''^ to AMM of the muon. We have demonstrated that this 
dependence is weak, in particular, in the region of energies up to s = 1.5 GeV^ the differences 
between predictions of ChPT and sQED are very small compared to the present experimental 
accuracy. In general, the deviation of ChPT from sQED increases with increasing the minimal 
photon energy Ecut- However this deviation becomes essential only if the energy of the photon 
exceeds 400 MeV, in the region where a'^f''^''^ itself is beyond the existing experimental precision. 
To observe such effects the experimental accuracy has to be increased by at least one order of 
magnitude. 

In fact, this small deviation is not surprising. Firstly, at fixed value of s the low-energy 
photon region is described similarly in both models and, as follows from the calculation, this 
region dominates in al^'^'"' . Secondly, the main contribution to the integral over s in Eq. pUll 
comes from the region of the p-resonance, which is accounted for in the same way in sQED 
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and ChPT through the VMD modeL Therefore, the integral oJ^'^'^ is not sensitive to the chiral 
dynamics (see also discussion in |21j). 

The developed approach has also allowed us to investigate the C-odd asymmetry of the cross 
section caused by the ISR - FSR interference. In general, measurements of the asymmetry can 
test the FSR amplitude. We considered radiation of the photon at the small angle relative 
to the direction of the electron momentum, ^ < = 15°. In these conditions the absolute 
value of asymmetry is of the order of 9q ~ 10~^. According to our calculations the difference 
between sQED and ChPT shows up only at the high photon energies uj > 0.3 GeV. For the 
smaller photon energies the difference is less than 1%. Since the asymmetry itself is less than 
10"^ for the selected collinear kinematics the experimental observation of this difference in the 
energy region u < 0.3 GeV is problematic. Thus the model dependence of the asymmetry can 
experimentally be observed only at close to the threshold region, 4m^ < 5'^ < 0.4 GeV^. 

To our view, the photon FSR from the two-pion channel in e~^e~ hadrons process shows 
the model dependence only near the two-pion threshold where the photon energy is large. In the 
bulk of energies (0.4 GeV^ < < s) scalar QED is sufficient to describe the FSR contribution 
to a^'^'"' and the C-odd asymmetry. In that way our results validate recent calculations 
performed in framework of sQED. 

It is plausible that the more complicated many-particle channels are more sensitive to the 
chiral dynamics. Another possibility to test chiral models is the region of the space-like photon 
momenta (Q^ < 0). In particular, the virtual Compton scattering on the pion, e~7r''^ — * e~7r^7, 
allows one to obtain information on the pion polarizabilities (see, e.g., [221 E]) ■ 
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Appendix A. General form of FSR tensor 

The amplitude of the reaction 7*((5) 7(/c)+7r+(p_|_)+7r^(p_) depends on three 4-momenta, 
which can be chosen as Q, k, and I = pj^ — p^. Here Q = pi +P2 is the total momentum of the 
e^e pair with Q"^ = s = AE'^. For on-mass-shell pions Q ■ I = k ■ I. In general, the second-rank 
Lorentz tensor M^'^{Q,k,l) can be decomposed through 10 independent tensors I24j : 

10 

M^'^Q, k,l) = Y, ^"HQ\ k\Q-k,k-l), (A.l) 

i=l 

= {^M^^ Q^^Q^^ ki^k", Pr, PQ^ Q^r, Fk", k^l", Q^'k^ k^Q"}, 

where parity conservation is taken into account. The tensor M^'^{Q,k,l) obeys the following 
properties: 

M^^iQ, k, I) = M'^'iQ, k, -I) = M^^i-k, -Q, I). (A.2) 

The first equality follows from the charge conjugation symmetry of the S-matrix element 
(7(fc), 7r"'"(p+)7r~(j9_)|5'|7*(Q)) = {'y{k),Tr~{p+)TT~^{p^)\S\'y*{Q)), and the second one is due to 
the photon crossing symmetry: Q —k and fi ^ u. Eqs. ()A.2|) impose certain constraints on 
the scalar functions Fi{Q'^, k'^,Q ■ k,k ■ I). 

The consideration below follows that of Ref. where the virtual Compton scattering 
7*(g) + 7c^{p) 7(5') + TT^{p') with the space-like initial photon (g^ < 0) and real final photon 
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(g'^ = 0) has been studied in detail. Some of the results for the reaction 7*(Q) —>■ 'j{k) + 

7r+(p+) +7r~(]9_) can be obtained from the corresponding results of [211 after the substitutions: 
pf^ ^ -pi^^^ p'M ^ ^ ^ QM^ q'f^ ^ ^m. 

The gauge invariance conditions Q^M^^{Q, k,l) = and M^'^{Q, k, l)ky = lead to the five 
linear equations between functions the Fi in Eq. flA.l|) . and in the general case of two virtual 
photons one is left with five scalar functions (see Eqs. (14) and (15) of [21] )• We are interested 
in the situation, where the final photon is real, i.e. A;^ = and k^e^ = [e^ is the polarization 
vector of the final photon), while the initial photon produced in e+e~ annihilation is virtual 
with > 4m^. In this case one obtains 

M^-^iQ, k, I) = -le^T^h + rr/2 + rt^'' h) ^ -ie'M^^{Q, k, I), (A.3) 

with the gauge invariant tensors 

T^'' = k-l{l^'Q'' -g^"'k-l)+F{k^k-l-l^k-Q), 

= Q\g^'''k-l~kn'')+Q^'{rk-Q-Q''k-l). (A.4) 

The scalar functions fi = fi{Q^, k ■ Q,k ■ I) are either even or odd with respect to the change 
of sign of the argument k ■ I: 

fi,2iQ^k-Q,k-l) = +fi,2iQ\k-Q,-k-l), 
h{Q\k-Q,k-l) = -MQ\k-Q,-k-l). (A.5) 

The factor — ze^ in Eq. ()A.3|1 is included for convenience. It thus follows that the evaluation of 
the FSR tensor amounts to the calculation of the scalar functions fi (i = 1, 2, 3). 



Appendix B. Chiral Lagrangian for pseudoscalar, vector 
and axial-vector mesons and photon 

We choose the 0{p^) chiral Lagrangian derived by Ecker et al. which contains vector 
mesons and axial-vector mesons. It includes interaction of pseudoscalar, vector and axial-vector 
mesons, and photon. The explicit form is given in ^1] : 

-^Tiiy^Vxf.V.V^ - \my^,Vn - iTr(V"AA^V.A'^^' - UnlA^^A^ 

+ ^Tr(V,./r) + '^T.{y,^u^u^) + ^Tr(A,./r), (B.l) 

where V = exp{i\/2^/ fT^), $ describes the SU{3) octet of pseudoscalar mesons, V^i, (A^iy) 
is antisymmetric field describing the SU{3) octet of polar-vector (axial-vector) mesons, and 
pfMu _ Qfij^u _ Qvj^fj, ^YiQ EM tensor, where the photon field is denoted by B^^. Further, 
Fy, Gy, Fa are constants, whose numerical values are specified in Sect. [H For more details 
on definitions and notation see Ref. [T^ . 

For treating the process 7* 7r+7r~7 at tree level it is sufficient to keep in Eq. ()B.1|) the 
terms containing the neutral meson p''(770), and the charged mesons a^(1260) and vr^, as well 
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as the photon. We obtain 



TT TT 



J TT 



(B.2) 



Lagrangian ()B.2j) leads to the Feynman rules discussed in Appendix C. The diagrams contribut- 
ing to the 7* 7r+7r~7 reaction at tree level are shown in Fig. 2. For a general case of the two 
virtual photons we obtain the FSR tensor Mp": 



1 



Q^-2Q- j9_ 



- 2F vG 

J TT 

2FyGy 



J TT 



T>'iQ~P-,P-)V'ip+,P+ + k) + (p+ ^ p_) - 2g>^'' 



1 1 

+ 



— Q'^ m?p — k"^ 



^"k-Q- k^'Q") + 



m?p — k^ 



^1 - (Q-p-) 



,[(k + p+rik ■ Q{k + p+r 



-k-{k + p+)) -Q-{k + p+){k^^{k + p+r -k-{k + p+)g>^n] + {P+ ^ P-] 
where the EM vertex for the off-mass-shell pion (with initial pi and final pj momenta) is 



(B.3) 



F G 

^ ^ [{Pi+PfYr^ -r^{p)-pl)l r = pf-pi. (B.4) 



fl{m? - P] 



Note that Lagrangian (jB.ljl was applied in Ref. in calculation of the Compton tensor for 
7*7r+ — > 7*7r+. 

Appendix C. Feynman rules for ChPT Lagrangian 

Following we describe the vector (axial-vector) meson by the antisymmetric tensor field 
that corresponds to the following form of the propagator 



[^"^/''(M^^ - q') + g^^qf^q" - g'^'q^q" - (a ^ /?)] 



(C.l) 



where M is a mass of the vector (axial-vector) meson. 

The vertices corresponding to the ChPT Lagrangian from Appendix B are listed in Fig. 7. 
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/ 

''''''' +(-)l^(fcy-^v^) 

5^ k •^^ '■h 



P- (f) 

-je(p+-p_)'' (a/3) + H^rs'" - 



^ a+(a-) ^ 



/ 
\ 



(f) 



(a/3) 



{a/3) '^'^7^ 



/ 

(,,) -v/^^'zzzzz (a/3) ^(/S(p_+p,)«-/"(p_+p,)3)-!|(r/''-fcV'') 



(a/j) = = = = = = (ri ^(/«(p_+p,)°-/>_+p,).«)+g(rs*-fr^9°'') 



Figure 7: 
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Appendix D. Expressions for functions /112 

The functions hi and /12 of Sect. 12.21 can be written as 
where h'f^^^ are determined from the following equations ^T] 



(D.l) 



(D.2) 



2hf^^ - hf'^^Q' = ■ Q? + - ml){-q' + (g^ - 2ml)L,)] (D.3) 



withe= 1 f) andLi = -ln- 1. 

V q^ J ^1-^ 

From Eqs. ()32|) - ()34j) we obtain the equation for A/i2 and (2A/ii — 5/^2) in ChPT 



2Ahi - sA/i2 = 2tt^ Ci + + 



k ■ Q [(m^ — m^) {r — k ■ Q) + 



Aml(k-QY- 



where 



{ml - ml){r - k ■ Q) 



Aml{ml + ml){k- QY , ^ 2 2 , / 2 2n/ , 0^2 , \ «^2 
[ml- m^)[r - k ■ Q) jk-Q 



Ci = 2(fc-Q)V' + a/(2g'fc-Q + rg2) + ia2((g2_2A;-Q)g2 + 2r2), 



(D.4) 



(D.5) 



Co 



- Y{(^ ■ + 4A; ■ Q + 4m2) + 2k ■ Q{Q^ + 2k ■ Q)r 

+ f{{k ■ Q)'[8ml + 2ml{Q' + 6k ■ Q) 

+k ■ Q{2Q^ + 5k ■ Q)] + 2k ■ Qr[-q^Q^ + 12k ■ Q{ml + ■ Q)] 
-2r\Q^ + mlq^ + k ■ Q[(f - \\k ■ Q)\ - 4QV=^ - 3r^}, 



(D.6) 



C3 = -{[8mt + 2m2(-g2 + 6A;-g) + 5(A;-g)2 + g4](^.g)2 + 2rA;-g 

X [g^ + 2Q^k ■ g + 4(fc ■ g)2 + 4m2 (A; ■ Q)\ + r2[g4 + sg^fc ■ g 

+2(A; ■ Qf + 2m^g'] + 4gV=^ + 



(D.7) 



(m^ — m;)(r — k ■ Q) 



18 



2,q\ml -ml){r-k-Q) + 2{k ■ Qfiq^ + 2ml) , 2_ 



SReFJs 



Ark ■ Q 

r{q^{ml - m^)(r - k ■ Q) + 4:{k ■ QYml) 



2- 



2k ■ Q{r — k ■ Q) {m\ — m, 
1 r + jk-Q 2 2 , n 



Fl r_ Fv- '^FyGy ( 1 1 
J X9 9 



LA, (D.8) 



Appendix E. Kinematics of the 3-particle final state 

Solving energy-momentum conservation for the pion energy in Eq. (fTBj) requires some care. 
Firstly, notice that the energy of the photon at fixed CM energy ^/s = 2E varies within the 
limits < uj < ujmax = E — ml/ E. Secondly, by requiring positiveness of the expression under 
the square root in Eq. (fT?)|) at arbitrary cos6'^+, we get the conditions 

2E{E-m^) 2E{E + m^) ,^ ^ , 

u < = , u > UJ+ = -. E.l 

2E-m^ - ^ 2E + m^ ^ ^ 

Clearly > uJmax and < uJmax- Thus the restriction for the photon energy is < u; < ct;_. 
The corresponding invariant mass of the tt'^7c~ pair is AE^ > g2 > 4:E^m^/{2E-m^) ^ 2Em^. 

The requirement < uj < uj^ coincides with the condition that the energy conservation 
low in Eq. ()14j) leads to one solution for E^, namely the one in Eq. ()16p . In the other case, if 
uj^ < UJ < Umax, the situatiou becomes more comphcated: there appear two solutions e\_ ^ and 
E^^ , which differ by the sign in front of the square root in Eq. ()16p . Correspondingly one has 
to sum over two terms in Eq. (fT^ corresponding to these solutions. Besides, the angle 6*^+ in 
this case is limited by the value 
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• 2. ^ _ AEiE-u)[EiE-u)-m^, 

Sm t^-y+)max — ^-^3 • {^■2) 



For these values of photon energies each angle 6^^ in the Lab frame (CM frame for colliding 
e"'"e~ beams) corresponds to the two different angles between momenta of tt^ and 7 in the tt+tt" 
CM frame. Here we refer to monograph [22] (ch. Ill), where these aspects of kinematics are 
considered in detail. 



Appendix F. Contribution to FSR from p tt 7 decays 

The diagrams with intermediate charged p-meson can be obtained from the 3rd row of diagrams 
in Fig. 2, if a^-meson is replaced by p^-meson. We choose the chiral Lagrangian, describing 
the odd-intrinsic-parity sector, in the form J7] 

L(^) = i^veM-/3Tr(Vn«^ /f }) ^ -^^^^e.^.^F^'' /Td-'Tf (F.l) 
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q^=0.4GeV'' 




40 80 120 160 

9^, deg 

Figure 8: Charge asymmetry as a function of pion polar angle for s = 1 GeV^. The solid 
(dotted) line corresponds to the tensor (vector) formulation for p-meson, the dashed line - 
calculation without p — > 7r7 contribution. 



in the vector formulation for the p-meson field, where e^j^a/j is the totally antisymmetric Levi- 
Civita tensor. 

The constant Hy can be determined from the p^ — 71^7 decay width r(p^ ir^'-f) = 68.7 
keV dni. From Eq. one finds 



27/2 ^ m2 



r(p^ - vr±7) = - (F.2) 



and thus Hy = 0.0363. 

The corresponding contribution to the invariant functions of Sect. Stakes the form 



8H, 



2 r 

V 



1 1 



C{1) C{-l)J \C{1) C{-1) 



(F.3) 



8Hl f 1 1 



^^^^ = ^(ciiy-cFo)' ^^-^^ 

where C{±1) = ml-{k + p±f -impTp{{k+p±f) with {k+p±f = {Q^ + 1^ + 2k-Q±Ak-l)/A. 

If we choose the antisymmetric-tensor field formulation for the p-meson as was done in the 
rest of this paper, then the Lagrangian, which is equivalent to ()F.1|) on the mass shell, reads 

L(-) = wF-^p-^vf. (F.6) 

For this Lagrangian the functions A/23 same as in ()F.4|) and ()F.5|) . while A/C differs 

from ()F.3|) by an additional term, 

iAfty = Aff + ^-^. (F.7) 

According to our calculations at invariant masses from the two-pion threshold to ~ 0.4 
GeV^ the p ^ 717 contribution to the charge asymmetry may be of the same order as the 
ai TT'-f contribution, if the tensor formulation for the p-meson field is applied (see Fig. 8). 
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For the higher values of the considered mechanism is suppressed with respect to other 
contributions. 

Regarding the seeming difference between vector and tensor formulations, we should note 
that, as argued in |T7l 12^], the effective Lagrangians in the two formulations would become 
equivalent if the Lagrangian in the tensor formulation included an additional local term. Ap- 
parently the contribution from this local term to the functions A/f would cancel the term 
64Hy/{9f^) in ()F.7|) making the charge asymmetry independent of the formulation for the 
p-meson field. Therefore we can conclude that the contribution of the 7* p^ir^ — > 7r"'"7r~7 
process to the asymmetry is very small at all two-pion invariant masses. 
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